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Abstract 



The Neel and the singlet RVB orders of the t-J model in a 2D square lattice 
are studied in the slave-boson mean-field approximation. It is shown that the 
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Neel order parameter takes the maximum value at the finite temperature and 
disappear at the lower temperature for a certain range of doping. It is also 
shown that the Neel and the singlet RVB orders coexist at low temperature. 
S This suggests the possibility of the coexistence of the Neel and the supercon- 

ducting orders. 
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U ■ I. INTRODUCTION 

In order to make clear the mechanism of high-T c superconductivity in cuprates, it is 
indispensable to understand the whole phase diagram on the temperature-doping plane. 
Recent theoretical studies by the slave-boson mean-field theory based on the t-J model show 
the validities of the method and the model for high-T c cuprates [ [1]— £|. These studies, 
however, do not take into account the Neel order explicitly. A few works based on the same 
theoretical scheme study the Neel order [ |5||J. In ref. the Neel order and the singlet 
RVB order parameters are considered, but the bond order parameters are not treated in 
the self-consistent way. In ref. |6] the Neel order and the bond order parameters of spinon 
and holon are included. It is shown that there is a phase transition between itinerant and 
localized antiferromagnetism in the Neel state. But the singlet RVB order is not studied 
there. 
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In this paper all of the order parameters, i.e. , the Neel order, the RVB order and the 
bond orders of spinons and holons are considered to study the whole phase diagram in the 
temperature-doping plane. It is found that the Neel and singlet RVB orders can coexist 
in the low temperature range. It is also found that the Neel order parameter takes the 
maximum value at the finite temperature and disappear at the lower temperature for a 
certain doping region. 

II. MODEL AND RESULTS 

In this communication we study the t-J model, which is given by 

H = ~t £ + c] a c ia ) + J £ Si ■ S j . (1) 

<i,j>,a <i,j> 

Here q ct = Ci a (l — ni- a ), where Cj CT is the electron annihilation operator at site i with spin a, 
Si is the spin operator at site i, t is the transfer integral, J is the exchange coupling constant 
and the summation with respect to < i,j > is taken over the nearest neighbor pairs. We 
employ the slave-boson method to solve Hamiltonian ([[]), in which Cj CT is replaced by the 
product of the spinon (fermion) and holon (boson) operators, f ia and bf c ia = fi a b\. The 
spinon and holon numbers should satisfy the local constraint, 

+ && = !■ (2) 

The form of the t-J model eq. ([!]) is different from that which is derived from the single- 
band Hubbard model in the large U limit. In the latter case the exchange term is JJ2{&i ' 
Sj — (l/4)rijnj). The both t-J model have a local SU(2) symmetry at half filling [ [TJ. If 
we approximate the local constraint (Q) to the global one, the latter t-J model might lose 
the symmetry. That is, however, preserved even in that approximation in the present case. 
There is another reason of omitting — (l/4)n.jrij in the J-term. In the case of hole doped 
high-T c cuprates it is considered that the original microscopic model is the d-p model. The 
t-J model is derived from the d-p model through the construction of the Zhang-Rice singlet 
[ [8|. In this case the magnitude of the coefficient in the —nirij term is different from 1/4 
and is estimated to be negligibly small [ |§. 

We treat the Hamiltonian in the mean-field approximation, where the local constraint 
is relaxed to the global one: J2i(J2a fio-fia + tybi) = N, where N is the number of the 
lattice sites. In order to consider the Neel order, we divide the square lattice into two 
sublattices A and B and define the Neel order parameter (staggered magnetization) m s = 
(SX) = —(Sb). We also take into account the uniform bond order parameters (J2 a fLfja) 
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and (b\bj), which represent the coherent motions of spinons and holons, respectively, and 
the singlet RVB order parameter {fiifi+(i,o)-\) and (fiifi+(o,i)-\), which represents the singlet 
pairing of spinons. We assume the ci-wave pairing: (/i|/i+(i,o)t) = ~~ (Al/«+(o,i)t)' shice the 
singlet RVB susceptibility for the rf-wave pairing is shown to diverge at the temperature 
higher than that for the s-wave pairing. 

The mean-field Hamiltonian is obtained as 

H M F — E ' \ ~~ V)fLfka + ( e fc+Q - ^fk+Qafk+Qa 
ka L 

- aJm s (f ka f k+ Q a + ft +Qa fka) 

./ 3 



+E' -A^f-nfU + A fe+Q /i fc _ Qi /i +QT + h.c.) 



+ E' 
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£ = UNA f A b + -JiVm 



- J7VA 2 + 3JNAI + AiV , 



(3) 

(4) 
(5) 

(6) 



and A is Lagrange multiplier, which corresponds a holon chemical potential, ji = /i/ + A 
is the spinon chemical potential,where \if is the electron chemical potential, Q = (7r, 7r), 
7| = 2(cos/c :c + cosfcj,), A fc = J(fiifi+(i,o)])li, it = 2(cos/c :E - cos^), and the prime 
on the summation symbol indicates that the sum is taken over the Brillouin zone for the 
sublattice. The diagonalized Hamiltonian has four energy bands, ±Ei k and ±i?2fc, where 
E lk = y/fe - /i) 2 + (3A fe /4)2 , E 2k = ^ fc + //) 2 + (3A fc /4)2 and ^ = + and 
/i < for 5 7^ 0. The self-consistent equations are given by 

1 1 
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E' 
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exp((co fe - A)/T) - 1 exp((u k+Q - A)/T) - 1 



1 E' 



tanh ■ 
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tanh 
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— E' 
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7fc+Q 



_exp((^ fc - A)/T) - 1 exp((w fc+ Q - A)/T) - 1 



— Y'lt- 
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(7-a) 
(7-b) 
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(Al/i+(l,0)T) = ^77T7 



32iV fc 



1 -^lfc 1 T^fc 

■ tanh— — + — — tanh- 



Ei k 2T E 2 k 2T 
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The above equations are solved numerically for t/J — 4.0, which is considered to be appro- 
priate for high T c cuprates. 

The resultant phase diagram is in Fig. 1. Here Td (dotted line), Tn (solid line) and 
Trvb (broken line) are the critical temperatures of the bond order parameters, the Neel 
order parameter and the singlet RVB order parameter, respectively, and in the meshed and 
hatched regions each order parameter is finite. Both of the two bond order parameters 
fiafjo) an d (b\bj) become finite at the same temperature Td above which the holons 
and spinons are localized and below which they are itinerant. In the left up hatched region 
enclosed by Tn line the Neel order exists. The RVB order survives below the Trvb line- 
In order to present the physical pictures for regions with different doping, we divide 
the T-8 plane into four regions of different 5. In the region of 5 smaller than 5± at which 
Tn equals Td, there are the two kind of antiferromagnetism below Tjsf. In the temperature 
Tn > T > Td, localized antiferromagnetism occurs which is due to localized spin. Below Td 
itinerant antiferromagnetism occurs. At Td there is the phase transition between localized 
and itinerant antiferromagnetism. At Trvb < Td the singlet RVB order coexists with the 
Neel order. It is to be noted that Trvb as well as Td go to zero for 5 = 0. In the region 
between 5\ and 8 C , spinons and holons become itinerant below Td and itinerant antiferro- 
magnetism occurs at Tn < Td- The singlet RVB order coexists with the Neel order below 
Trvb < Tn. The region between 8 C and 62 is quite interesting. In this region the Neel order 
parameter exists only in the finite temperature range. The physical reason for this behavior 
will be discussed in the next section. At higher doping rate 5 > 62, Trvb decreases mono- 
tonically with 5. At 5 > 0.275 Trvb is lower than 10~ 3 . In principle, Trvb is finite for finite 
spinon density. We could not determine the finite value of Trvb for the sake of numerical 
accuracy. 

The temperature dependence of the staggered magnetization m s is shown in Fig. 2 for 
various doping rates. The order parameter takes the maximum value, i.e. , the hump at 
< T < Tn . This hump is particularly noteworthy for high doping. For the doping rate 
5 C < 5 < 82, the order parameter, which increases once below Tn, vanishes at the finite 
temperature. The Neel order parameter for a fixed temperature decreases as the doping 
increases. 
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III. DISCUSSION AND SUMMARY 



We have studied the t- J model on the 2D square lattice in the mean-field approximation 
with the slave-boson description. We have taken account of four order parameters, viz. 
, the Neel order parameter, the d-wave singlet RVB order parameter and the bond order 
parameters of spinons and holons, and have obtained the phase diagram in the T-S plane. 

Our result of the temperature dependence of the staggered magnetization m s in Fig. 2 



1C] agrees qualitatively with the neutron scattering experiments of L^CuO^ by Keimer 



et. al. [ pn[| . The staggered magnetization curve against temperature has pronounced peak 
for samples with lower Neel temperature [ In their experiments the doping rate of each 
samples (Fig. 2(a) in their paper) is not clear, but the sample with lower Neel temperature 
is considered to be the sample with the larger doping rate. Therefore their experimental 
findings is consistent with our result Fig. 2, where the peak of the staggered magnetization 
is pronounced for larger doping rate. 

The physical reason why the staggered magnetization is reduced at the low temperature 
side is the following. In general the nesting of Fermi surface with fixed wave vector Q = (n, n) 
is hindered by the doping, since the Fermi surface is shifted from that at half filling, and 
susceptibility for the Q is reduced. This reduction is pronounced at lower temperature since 
at the higher temperature the effect of the shift of the Fermi surface from that at half filling 
is smeared out [ 0,f| • 

In Fig. 1 the Neel order coexists with the singlet RVB order for the low temperature 
side. This suggests the coexistence of the the Neel order and the superconducting order. 
The superconducting order is characterized by the coexistence of the singlet RVB state and 
the bose condensation of holons because of (cjjCjf) = (fiifj\)(bibj). The latter does not 
occur in the exact 2D system at finite temperature, but may occur if 3D effect is taken into 
account in the holon dispersion. Then the superconducting order can coexist with the Neel 
order in the coexistent phase of the Neel and singlet RVB orders. This coexistent phase has 
not been observed in the experiment. The reason is considered that the localization effect 
induced by disorder associated with doping is much enhanced in the small doping region 
in the Neel order phase, because the effective bandwidth and the carrier number are quite 
small there. We expect the experimental observation of the coexistence phase in samples 
where the disorder effect is small. For example La2Cu04 +( 5 may be a good candidate since 
the metallic state with the Neel order is observed and the effect of disorder may be weak. 
The preparation of the sample with single phase are desired. 

In summary, we have studied the 2D t- J model in the slave-boson mean-field approxima- 
tion with four order parameters: the Neel, RVB, and bond orders and the phase diagram has 
been obtained. It has been shown that the Neel and the singlet RVB orders can coexist at 
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low temperature. We also show that the staggered magnetization takes the maximum value 
at the finite temperature. The interesting reentrant behavior of the staggered magnetization 
is predicted for a certain range of 5. 
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FIGURE CAPTIONS 



Fig. 1. The ^-dependence of the critical temperatures for t/J—A.O. TD(dotted line), 
Trvb (broken line) and (solid line) are the critical temperatures of the bond order of 
spinon and holon, the RVB order and Neel order parameter respectively, and in the meshed 
and hatched regions each order parameter becomes non-zero. In the overlapping regions, the 
order parameters coexist. At Si (~ 0.02) T D equals T N . At 5 C (~ 0.15) T N =0. 5 2 (^ 0.16) is 
the maximum doping rate where the Neel order exists. 

Fig. 2. The temperature dependence of the Neel order parameter, i?in, at various doping 
rates. 
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